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Abstract . We have used two years of strain 
tide data to study the response of the Earth to 
the diurnal and semidiurnal tidal excitations. 

Our results show that there is significant struc- 
ture in the response of the earth to tidal exci- 
tations near one cycle/sidereal day. This struc- 
ture agrees with the resonance behavior predicted 
from the calculations of the forced elastic- 
gravitational response of an elliptical, rotating 
earth with a liquid outer core. The data can 
also be used to test for possible preferred 
frames and spatial anisotropies. We find that 
upper bounds on the parameterized post-Newtonian 
(PPN) parameters which characterize these effects 
are a-2 - 0,007 and S 0.005. 

Introduction 

We have analyzed approximately two years of 
strain-tide data obtained using the 30-meter 
laser strainmeter we have previously described 
(Levine and Hall, 1972). 

Data Acquisition 

The data were obtained using a 30-meter long 
laser strainmeter located in the Poorman Mine, 
an unworked gold mine located approximately 8 km 
west of Boulder, Colorado at latitude 40.03“N, 
and longitude 254.67®E. 

The heart of the strainmeter is an evacuated 
30-meter Fabry-Perot Interferometer located along 
the length of the tunnel. The axis of the inter- 
ferometer is 7° west of North. 

The interferometer is illuminated bya3.39-um 
helium-neon laser. A servo loop piezoelectri- 
cally tunes the laser to keep its wavelength 
coincident with one of the transmission maxima 
of the long interferometer. The frequency of 
the laser, f , is therefore related to the length 
of the interferometer, L, liy 



where n is an Integer and c is the velocity of 
light . Thus 
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A second 3.39-pm laser is stabilized using 
saturated absorption in methane. The beat fre- 
quency between the two lasers is extracted for 
further processing. Then 
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so that a measurement of the fluctuations in the 
beat frequency provides a direct measurement of 
the fractional change in the length of the long 
path. The relationship between beat frequency 
and strain has no adjustable constants or cali- 
bration factors. 

The beat frequency is digitally recorded 10 
times/hour along with other information including 
local barometric pressure, etc. The digital data 
are bandpass filtered using a symmetric convolu- 
tion filter and then decimated to one sample 
every two hours for comparison with theory. 

Data Analysis 

We have used all of the components published 
by Cartwright and Edden (1973) in our analysis. 
For each component with frequency fj^ we construct 
a time series of the form 


a^ cos (2tfj^t +(j>j^ + aj^)Y™(0,<p)T(n,m,e,<C,0g) 

where aj^ is the amplitude given by Cartwright and 
Edden, ()>ij is the phase of the component, Y® is 
the spherical harmonic computed at the station 
co-latitude 0 and East Longitude q>. The quantity 
T is a function converting potential to strain 
along the axis of the strainmeter 0g (Levine and 
Harrison, 1976). The quantity is -90° if (n+m) 

is odd and is zero otherwise (Cartwright and 
Tayler, 1971). 

The theoretical series were further modified 
by a function to correct for local topography, 
local crustal inhomogeneities, cavity effects and 
ocean loads (Levine and Harrison, 1976). We as- 
sumed initially that these effects do not vary 
rapidly with frequency, so that all of the diur- 
nal components have the same correction as Oj , 
and that all of the semidiurnal components have 
the same correction as M 2 . 

The terms are then grouped by frequency. Each 
frequency group contains all of the terms (re- 
gardless of parentage) which differ by less than 
one cycle/yeaf from each other. 

This process produces 48 time series. We fit 
these series to our data by the least-squares 
method using an adjustable amplitude and an ad- 
justable phase for each cycle/year group. 

The results of this process are shown in Fig. 1 
for the diurnal amplitude. Fig. 2 for the diurnal 
phase, and Fig. 3 for the semidiurnal amplitude 
and phase. The error bars represent one standard 
devi ation and are obtained from estimates of the 
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Fig. 1. Normalized transfer function for the 
diurnal tides. For clarity only the major com- 
ponents are plotted. ' The error bars are one 
standard deviation. 
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Fig. 2. Normalized transfer function phase for 
the diurnal tides. For clarity only the major 
components are plotted. The error bars are one 
standard deviation. 
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Fig. 3. Normalized transfer function amplitude 
(lower curve) and phase (upper curve) for the 
semidiurnal tides. For clarity only the major 
components are platted. The error bars are one 
standard deviation. 
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Fig. 4. Normalized transfer function amplitude 
for the diurnal tides when ' the frequency-depen- 
dent Love numbers are Inserted into the fitting 
function. The dashed line is obtained by fitting 
a straight line to the transfer function ampli- 
tudes weighted by their respective uncertainties. 



Conclusions 


noise. 

The amplitude of the diurnal transfer function 
shows a statistically significant structure. 

There appears to be a significant decrease in 
the transfer function for the lower frequency 
diurnal components. We attribute this to a 
slow change in the contribution of the ocean 
load to our observed data. 

More significantly there is a dip in the 
transfer function near one cycle/day (we have 
not plotted the Sj amplitude at one cycle/day 
since it is heavily contaminated by thermoelastic 
processes). This dip is consistent with the ef- 
fects of the nearly diurnal resonance associated 
with the liquid core. The consistency can be 
shown more clearly by inserting the resonance 
directly into the fitting function. This may be 
done using the frequency dependent Love numbers 
published by Molodensky (1961) and Shen and 
Mansinha (1976). These Love numbers produce 
changes in the transfer function at all of the 
diurnal frequencies but, except near one cycle/ 
day, the new transfer function lies within one 
standard deviation of the old one. The modified 
transfer function is shown in Fig. 4. As can be 
seen from the figure the resonance models account 
for the dip in the transfer function near one 
cycle/day. Unfortunately, the data cannot be 
used to infer the fine structure near the reso- 
nance. Furthermore, the Insertions of realistic 
estimates for the energy loss in the nearly diur- 
nal band may affect the shapes predicted by all 
of the models somewhat. It Is unlikely we will 
be able to test any of these effects using our 
data. The strongest test for the resonance ob- 
tainable from our data comes from the components 
Pi and Ki which are really on the tail of the 
resonance function. Our measurements at 4n and 
4>l (where the resonance has the largest effect) 
have error bars whose size probably does not per- 
mit them to be used in a quantitative comparison 
between theory and experiment or in attempts to 
differentiate among the various calculations of 
the effect of the resonance. 

These results may be compared with the analy- 
sis of Warburton and Goodkind (1977) using data 
from a cryogenic gravimeter, with the gravimeter 
data analysis of Abours and Lecolazet (1978) and 
with the analyses of Lecolazet and Melchior (1977). 
Their results are generally in agreement with 
ours, confirming the general shape of the reso- 
nance. None of the analyses is able to make a 
quantitative comparison with theory because of 
the relatively poor signal-to-noise ratio in the 
measurements of the crucial components <tii and i|<i. 

The semidiurnal transfer function amplitude 
shows far less structure, and the agreement be- 
tween experiment and theory is quite good. There 
is no evidence of anomaly at $2 (two cycles/day) 
confirming that the anomaly at Sj is almost cer- 
tainly of thermoelastic origin. 

We may place upper limits on the various PPN 
parameters by calculating the magnitudes of the 
anomalous tidal components in terms of the PPN 
parameters 02 and In this way we conclude 

02 S 0.007 and C s 0.005. 


These results confirm the general correctness 
of the various published earth models. From the 
point of view of the current discussion they are 
also significant in that they illustrate the 
sort of measurements that can be performed with 
laser strainmeters. These measurements show 
that we can make meaningful measurements of 
diurnal strain changes at the 10 level, and 
that almost continuous operation for several 
years is possible. 

From the geodetic point of view, the limitation 
on the utility of laser strainmeters arises from 
their sensitivity to local effects, especially to 
spurious motions near the piers . 

Spurious motions of the piers at the level of 
millimeters will even play a significant role in 
geodetic measurements made over much longer base- 
lines. Motions of this magnitude represent frac- 
tional changes of parts in 10® even over 50 km 
baselines, so that such effects will make sig- 
nificant contributions to the error budget of 
any geodetic instrument now in operation or under 
construction. 

It is important to compare the data obtained 
with laser strainmeters with measurements obtained 
using electromagnetic distance measuring equipment 
operating over parallel baselines. If this com- 
parison shows that the two provide a consistent 
picture of the regional strain field, then laser 
strainmeters may prove useful in measuring re- 
gional strain at fixed observatories. In this 
service they can provide significantly higher 
sensitivity than any other technique. 
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